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AN IMPROVED PURE SOURCE TRANSFER DOMAIN 
DECOMPOSITION METHOD FOR HELMHOLTZ EQUATIONS IN 

UNBOUNDED DOMAIN 

YU DU* AND HAIJUN WU* 

Abstract. We propose an improved pure source transfer domain decomposition method (PST- 
DDM) for solving the truncated perfectly matched layer (PML) approximation in bounded domain 
of Helmholtz scattering problem. The method is based on the the source transfer domain decompo¬ 
sition method (STDDM) proposed by Chen and Xiang. The idea is decomposing the domain into 
non-overlapping layers and transferring the sources equivalently layer by layer so that the solution 
in each layer can be solved using a PML method defined locally outside its two adjacent layers. 
Furthermore, we divide the domain into non-overlapping blocks and solve the solution in each block 
by using a PML method defined locally outside its adjacent blocks. The convergence analysis of the 
method is provided for the case of constant wave number. Finally, numerical examples are provided 
where the method is used as both a direct solver and an efficient preconditioner in the GMRES 
method for solving the Helmholtz equation. 

Key words. Helmholtz equation, large wave number, perfectly matched layer, source transfer, 
domain decomposition method 
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1. Introduction. This paper is devoted to a domain decomposition method 
based on the STDDM method (cf. [TK] i for the Helmholtz problem in the full space 
R 2 with Sommerfeld radiation condition: 

(1.1) A u + k 2 u = f in R 2 , 

= o(r~ 1 / 2 ) as r = |ar| —» oo. 

where the wave number k is positive and / G R 1 (R 2 ) 7 having compact support, where 
H 1 ( R 2 )' is the dual space of iU(R 2 ). The problem is satisfied in a weak sense (cf. 
[23]). 

Helmholtz boundary value problems appear in various applications, for example, 
in the context of inverse and scattering problems. Since the huge number of degrees of 
freedom is required resulting from the pollution error and the highly indefinite nature 
of Helmholtz problem with large number wave [3 [21131 EH 1IE1 [HI HH [Ml [Ml (Ml (Ml 
[251 [51] , it is challenging to solve the algebraic linear equations resulting from the finite 
difference or finite element method with large wave number. Considerable efforts in 
the literature have been made. One way is to find efficient and cheap methods [a m 
HZJEUE21I22, such as the coutinous interior penalty finite element method ['). 36, 36! 

, which use less degrees of freedom as the same relative error reached. Another way 
is to find efficient algorithms for solving discrete Helmholtz equations, e.g. Benamou 
and Despres [2], Gander et al [23] for domain decomposition techniques and Brandt 
and Livshit 0, Elman et at m for multigrid methods. Recently Engquist and Ying 
constructed a new sweeping preconditioner for the interior solution [201 . Then Chen 
and Xiang proposed the source transfer domain decomposition method (STDDM) 
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HE in which only some local PML problems defined locally outside the union of two 
layers are solved. Thus the complexity of STDDM is the sum of the complexity of 
the algorithms for solving those local problems which reduce the complexity to solve 
the whole linear system. We are inspired by the key idea of STDDM, and the main 
lemmas and idea of proofs also come from their work. 

In this paper we show the improved pure source transfer domain decomposition 
method (PSTDDM). 

Let Bi = {x = ( Xi,x 2 ) t £ R 2 : \x\\ < h, \x2\ < h}- Assume that / is supported 
in Bi. We divide the interval (—I 2 J 2 ) into N segments with the points (i = —I 2 + 
(i — 1)A£ where = 2 I 2 /N. Then we denote the layers by 

D 0 = {x = (xi,x 2 ) T £ R 2 : x 2 < Ci}, 

Dj = {x = (xi,x 2 ) T £ R 2 : C i <x 2 < Ci+i}, i = 1, - • • , N, 
il N +i = {x = (xi,x 2 ) T £ R 2 : Cjv +1 < x 2 }- 

Clearly, supp / C U^ 1 D i . Let fi = f in fl, and fi=0 elsewhere. Let = fi and 
In = /jv- The key idea is that by defining the source transfer function T* in the 
sense that 

[ fHy) G ix,y)dy = [ y+ +1 (f+)(y)G{x,y)dy Vx£U^++ 2 D,, 

J J Q i+1 

[ Ii{y)G{x,y)dy= [ ^r_ 1 (fr)( y )G(x,y)dy 

J Qi 

and letting f^_ 1 = f i±1 + then we have for any x £ Dj 

(1.3) u(x)=f-f fi(y)G{x,y)dy — f f+_ 1 (y)G(x,y)dy\ + 

\ J J£ii— 1 / 

(- [ fr+i(y) G ( x >y) d y\ 

\ Jn i+1 j 

Here, G{x, y) is the Green’s function of the problem fLT|)-([L2|. Observing 
know that u(x) in fli consists of two independent parts. The first part only involves 
the sources in fli and fb_i and the second one only involves the source in O l+ i. Thus 
they could be solved independently by using the PML method outside only Dj_i UDj 
and fli U fb + i respectively. Similar to STDDM, our PSTDDM also consists of two 
steps which could run in parallel and the complexity of every step is the same as 
that of STDDM. The error of PSTDDM are not larger than that of STDDM if the 
same numerical algorithm, such as the finite element or difference method, was used. 
Besides, since every step of PSTDDM just consists of local PML problems defined 
in the whole space, not in the half-space with Dirichlet boundary, these local PML 
problems also can be solved by using our PSTDDM recursively. As a result, the 
computational domain is divided into blocks and the local PML problems needed to 
be solved are defined outside the union of four adjacent blocks. 

The perfectly matched layer (PML) is a mesh termination technique of effective¬ 
ness, simplicity and flexibility in computational wave propagation. After the pioneer¬ 
ing work of Berenger mm, various constructions of PML absorbing layers have been 
proposed and many theoretical results about Helmholtz problem, such as those about 
the convergence and stability, have been studied pj m m ceei uni na up. In this 
paper, the uniaxial PML methods will be used. 


2.1), we 
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The remainder of this paper is organized as follows. In section [2j pure source 
transfer domain decomposition method in R 2 and some important lemmas and theo¬ 
rems, which are fundamental and illuminating for the PSTDDM in truncated domain, 
are introduced. Section [T3] shows the algorithm in the truncated bounded domain and 
the main result, that is, the exponentially convergence of the solution of PSTDDM 
in the truncated domain to the solution in R 2 . In section [ 3 J we show that the local 
PML problems can be solved by out PSTDDM recursively and proof the convergence 
of the method. At last, we both use the method as a direct solver and an efficient 
preconditioner in the GMRES method to solve the Helmholtz equation. 


2. Source transfer layer by layer. In this section, we introduce the PSTDDM 
for the PML method in the whole space. First, we recall the progress of deriving the 
PML method and set the medium properties of perfect matched lays which are a 
bit different from traditional medium and would be used in the following lemmas and 
theorems m■ In subsection tOl we also recall some basic lemmas. Then we introduce 


the two steps of the pore source transfer domain decomposition method in 


2.1. The PML method. In this subsection, we recall some knowledge about 
the PML method. 

The exact solution of equation ([IT]) with the radiation condition 1.2 can be written 
as the acoustic volume potential. Let G[x , y) be the fundamental solution of the 
Helmholtz problem 


AG(x, y ) + k 2 G(x, y) = —S y (x) in 


We know G(x, y) = \H^\k \x — y\) where H^\z), for z & C, is the first kind Hankel 
function of order zero. 


Then, the solution of (1.1) is given by 


( 2 . 1 ) 


l (x) = - f f(v)G(x,y)dy Wx 
Js . 2 


In this paper we used the uniaxial PML method 0113 mm] . the model medium 
properties are defined by 

ai(ar) = 1 + i<Ti(a:i), a 2 (x 2 ) = 1 + icr 2 (x 2 ) 

(Jj{t) = crj(—t) for t £ R 2 , <jj = 0 for |f| < lj, Uj = y 0 > 0 for \t \ > lj. 

where aj{xj) £ C 1 (R 2 ) are piecewise smooth functions, lj > lj is fixed and 70 is a 
constant. 

For x = (xi,x 2 ) T , we define the complex coordinate as x(x) = (ii(a;i),a? 2 (:ri)), 
where 

r x o r x o 

(2.2) Xj(xj) = / aj(t)dt = Xj+ i / (Tj(t)dt , j = 1,2. 

Jo Jo 

We remark that this kind of definition has been proposed in [15] and recall that the 
requirement, 07 = 70 for |f| > lj , is very important because of the use of proving the 
local inf-sup condition ( |2.33[ ) (cf. |T5]) for the truncated PML problem by using the 
reflection argument of mm and estimating the dependence of the inf-sup constants 
on the wave number k. 
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The complex distance is defined as 

(2.3) p(x,y) = [(iiOri) - yi{yi)) 2 + {x 2 {x 2 ) — ^(ya)) 2 ] 1 “ ■ 

Here, 2 1 / 2 denote the analytic branch of yfz such that Re( 2 1 / 2 ) > 0 for z £ C\[0, + 00 ). 
The solution to the PML problem is 


(2.4) u(x) = u(x) = - f f{y)G(x,y)dy Vx 6 R 2 , 

J R 2 

Since / is supported inside Bi , we know that y = y and u = u in Bi. 
The solution (2.4) satisfies the PML equation 

(2.5) J _1 V • (AVu) + k 2 u = f in R 2 , 


which could be obtained by the fact that Ats + k 2 u = / in R 2 and using the chain 
rule, where A(x) = diag and J ( x ) = oi 1 {x 1 )a 2 {x 2 )■ 

The weak formulation of ( |2.5[ ) is: Find u £ H 1 (R 2 ) such that 

(2.6) (AV-u, Vv) — k 2 (Ju, v) = — (Jf, v) \/v £ H\ K 2 ), 


where (•,•) is the inner product in L 2 (R 2 ) and (•,•) is the duality pairing between 
H 1 (R 2 )' and H\R 2 ). 

We have the following inf-sup condition for the sesquilinear form associated with 
the PML problem in R 2 which has been proved (cf. jT5j, Lemma 3.3): 


(2.7) 


sup 


I(HV<?!>, VVO -k 2 {J(j)^)\ 


ll^ll 


HH 


> Mo ||</>|| ff i( K 2) V(t>£H 1 { R 2 ), 


where the inf-sup condition p 0 1 < Ck 3 ^ 2 which is fundamental to our estimates. 
The fundamental solution of the PML equation (2.5) is (cf. [71 131|! 


( 2 . 8 ) 


G(x,y) = J(y)G(x, y) = ^ J(y)H i Q 1 \kp(x,y )). 


2.2. The PSTDDM for the PML equation in R 2 . In this subsection, we 
introduce our PSTDDM for the PML equation in the whole space and give the fun¬ 
damental theorems. 

We recall that the domain {.t = (xi,x 2 ) : \x 2 \ < Z 2 } is divided into N layers (??) 
and that fi{x) = f{x)\n i for any x £ fli and fi(x) = 0 for any x £ R 2 \fV We define 
smooth functions ^{x 2 ) and /3~ (x 2 ) by 


(2.9) 


Pt = !, Pi = 0, Pf = Pi ' = 0 as x 2 < Ci, 
Pt = 0, P~ = 1, Pf = Pf = 0 as x 2 > Ci+i, 



< C'(vc)- 1 , 



< C'(vc)- 1 , 


where C is a constant independent of Q, Q +1 and the subscript i. Our PSTDDM 
consists of two steps |Tj and [2j Clearly, the two steps [T] and [2] are independent of each 
other and can be computed in parallel. 











pure source transfer domain decomposition method 


5 


Algorithm 1 Source Transfer I for PML problem in R 2 

1. Let ft = h\ 

2. While i = !,••■ ,N - 2 do 

• Find uf € H 1 {R 2 ) such that 

(2.10) J~ X V ■ (AVu+) + k 2 u+ = ~ft~ f i+1 inR 2 

• Compute 

(2.11) ^ + +1 (/+) = J-'V • (AV(p++ lO) + k 2 Wt+iO- 

• Set 

( 2 - 12 ) f&i = fi+i+*t+iV?) 

in Oj + i and ft 1 = 0 elsewhere. 

End while 

3. For i = N — 1, find u't_ 1 € iJ 1 (M 2 ) such that 

(2.13) J- 1 ^ ■(AVu+_ 1 ) + k 2 u+_ 1 = -f+_ 1 -f N in R 2 


Algorithm 2 Source Transfer II for PML problem in R 2 

1. Let /jy = /tv; 

2. While i = N, ■ ■ ■ ,3, 

• Find u~ € U 1 (R 2 ) such that 

(2.14) J^V ■ ( AVu~) + k 2 u~ = -fr in R 2 

• Compute 


*T-i(fr) = J_lv • (AV(j9r lU r)) + k 2 0 Sr lU r). 

• Set f~_ x = /j_i + 'F“- 1 (/r) in f^i-i and ft_ x = 0 elsewhere. 
End while 

3. For i=2, find ulf € ff 1 (R 2 ) such that 

(2.15) J _1 V • (AV 1 I 2 ) + k 2 U 2 = -/1 - / 2 " in R 2 . 


By (2.10), (2.14) and (2.15), we know that uf is given by 

{ft + fi+i)J(y)G(x,y)dy \/x € R 2 , * = !,•■■ ,N- 1, 




(2.16) uf(x) = 
and u~ is given by 

(2.17) u~{ x)= [ f~ {y)J{y)G{x, y)dy \/x G R 2 , i = ,3 


(2.18) 


u 2 0*0 = / if 2 (y) + fi{y))J{y)G{x,y)dy. 


> UQ2 
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By simple calculation, we have the equivalent form of the source transfer operator 

\|/+ • 

W i+1- 

(2.19) ^ + i(A + ) = J- 1 ^ (AVf3t +1 u+) + ■ {AVuf) - Pt+Ji+i- 

and it’s easily obtained that 4+ /^i+i/i+i is in L 2 (fl i+ i) and supported in 
f2i+i. Similarly, we can get the equivalent form for \E')d 1 : 

(2.20) = J " lv iAVPZ. i«r) + J- x V^r , • (AVuT). 

The proof of the following two lemmas is quit similar to Lemma 2.6 in |15j . We 
omit the details. 

Lemma 2.1. Fori = !,■■■ ,1V —2, we have uY £ H 1 (M 2 ) and || M i~ ||/fi(R 2 ) — 
C ■ Let Mo = l, Mi = y/2M + (1 + y/2)Mi_i, where l is the diameter of Bi 

and M = max(Zi,/ 2 )- Then there exists a constant C > 0 such that 

k+(z)| + |Aw+(x)| < Ce~^ kl0 ^ ||/[| H i (Bi) , V]ar| > M u 

where is the dual space of H 1 (Bi). 

Lemma 2.2. For i = N, ••• ,3, we have u~ £ B 1 (R 2 ) and ||ti” H^i^) < 

C \\f\\ H i(Bi)' ■ Let Mo = l, Mi = \/2 M + (1 + \/2)Mi-\, where l is the diameter 
of Bi and M = max(Zi,l 2 )- Then there exists a constant C > 0 such that 

\uf(x)\ + |Au“(x)| < Ce“s fc7o|x| \\f\\ H i {Bl y V|x| > M u 


where Ft 1 ( Bi)' i s the dual s pace of Fl 1 (Bi). 


Lemma 


2.1 


and Lemma 


2.2 


show that uj and iq decay exponentially at infinity, 


which will be used in the following theorems. 

Theorem 2.3. The following assertions hold: 

(i) For i = 1, ■ ■ • , N — 2, we have, for any x £ fl((Ci+ 2 , +oo), 


( 2 . 21 ) 


ft(y)G{x,y)dy = / ^+(f+){y)G(x,y)dy. 


(ii) For the solution uf in (2.10), we have, for any x £ i = 1, • ■ • , N — 1, 

( 2 . 22 ) 


u t{ x )=[ f{y)G{x,y)dy. 

J fi(-00,Ci + 2) 


Proof. We first prove (|2.21|). By the property of (|2.8|) (cf. e.g. [[15], 2.11-2.13] 

[ESI T 

and y £ fl(Ci,Ci+ 2 ) 


[T, Theorem 2.8] and [[31]. Theorem 4.1]), we know that for any x £ fl(C,; + 2 ,+oo) 


v„ • {AVy( J-'Gix, y))) + k 2 J( J-^ix, y)) = 0. 

For x £ fl(£i+ 2 , +oo), y £ VLj, j = 1, • • • ,i + 1, G{x, y) decays exponentially as \y\ —» 0 
(cf. [15], Lemma 2.5). By Lemma 2.1 we know that uf(y) decays also exponentially 
at infinity. By integrating by parts, we have 

[ f?G{x,y)dy 

Jfi, 

= ~ [ ■ {AV v uf{y)) + k 2 Ju+(y)]G(x,y)dy 

Jfi(-oo,Ci+i) 


L 


(AVyUf {y) ■ e 2 )J G{x, y) - (AV y (J G{x , y)) ■ e 2 )uj ( y ) ds{y), 
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where ei is the unit vector in the X 2 axis. By using (2.91, we can do integration by 
parts to have 

[ f?G{x,y)dy= f [(AV y (^ +1 uf(y)) ■ n) J _1 G(:r, y)- 

J fli J dfli- 1-1 

(AV y ( J-* l G(x, y)) ■ n)Pt +1 uj (y)] ds(y) 

= [ J- l [\7 y • (. AV y (f3+ +1 ut(y ))) + k 2 Jpt + iul(y)}G(x,y)dy 

*/Qi+i 

= / ' K+i(ft)(y)G(x,y)dy, 

where n is the unit outer normal to 

Since y(y) = y and J{y) = 1 for any y £ Bi, By using (2.161 and (2.21) we could 
prove (2.22). For any x £ 0,; + i 


u ti x )= (ft + fi+i)J(y)G(x,y)dy 

= [ fi+i(y)G(x,y)dy+ [ fi{y)G(x,y)dy + [ ^f(f^L 1 )(y)J(y)G(x,y)dy 




f(y)G(x, y)dy 


' 1 


ft-i(y) J (y)G(x,y)dy 


f(y)G(x,y)dy. 


J u -ii^ 

This completes the proof. □ 

The second step [2] is similar to the first one |T] of the PSTDDM for the PML 
equation in the whole space. So by argument similar to the proof above, we can easily 
obtain the following results. 

Theorem 2.4. The following assertions hold: 

(i) For i = N,- ■ ■ ,3, we have, for any x £ $7(— 00 , Ci-i)> 


(2.23) 


fi (y)G(x, y)dy = / (/; )(y)G{x,y)dy. 

ja_. 


(ii) For the solution u i , i = N, ■ ■ ■ ,3, in (2.14), we have, for any x £ n,_i, 


(2.24) 


h 0 ) = / 

Jn 


n(Ci,+oo) 


f(y)G(x,y)dy. 


(iii) For the solution u 2 (x) in (2.15), we have, for any x £ fii, 

(2.25) u 2 {x)= ( f(y)G(x,y)dy. 

J R 2 


Combining Theorem |2.3| and Theorem |2.4| we could obtain the main result in this 
section. 

Theorem 2.5. We define u(j~(x) = 0 and ujj +1 (x) = 0 for any x £ K 2 . For any 
x £ f li, i = 1, • • • ,N, we have 

u{x) = -(uf^x) +u~ +l (x)). 
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Proof. From 
definition of u(x) 


(2.25), it’s easy to see that the lemma holds for i = 1. Using the 
(2.4) and (2.22), (2.24), we have, for any x G flj, i = 2, • • • , iV, 


u(x) = ~ f f(y)G(x, y)dy 

Jr 2 


= -[ f{y)G{x,y)dy- 

\JCl(-OoXi + l) 

= —(,uf-i(x) + uf +1 {x)), 


/fi(C;+i,+oo) 


f(y)G(x, y)dy 


where we have used y(y) = y in Bi. □ 

2.3. The PSTDDM for the PML equation in the truncated bounded 
domain. The PSTDDM for PML equation in the truncated bounded domain and the 
most important results in this paper are introduced in this section. First we introduce 
some notation. Let U be a bounded domain in R 2 and dU = T. Then the weighted 
norms are written as 


IMIz/qc/) — (1I^ 7 u IIl 2 (!7) + \\^ u \\l 2 (U) 
where du = diam(U) and 


1/2 


Mlj? l /*( r ) — Me/ 


Ml W) + Mf 


' r ) 


1/2 


I Ur 


I’d(cc) — v(x')\“ 


ds(x)ds(x'). 


/r Jr \x — x' | 

The following inequality are given (cf. [[T5j, 3.1]), 

(2.26) |MMi/ 2 (r) < (Ml du 1 ) 1 / 2 |M| L „ (r) + |F| II Vw|| L . (r) Vw € LF li00 (r),, 


The inequality (2.26) is easily derived from the definition of weighted norms. 
F< 

(2.27) 


For simplicity, the following assumption about the medium property is adopted: 

^i+<^2 /*^2 ~\~d>2 


ri 2-h«2 /'(-l-tai 

HI / cri(t)dt = / <j2(t)dt =: a, / o\(t)dt > d 

dll ^2 d ll~\~d2 

and 1 1 < I 2 , d\ = 2 c? 2 . 


This assumption is not essential. Those lemmas and theorems are also valid with 
a bit modification of the proof if the assumption is changed. 

Denote by Bl = (~h — di,U + d\) x (—Z 2 — d2,h + ^2) where l\ + d± > l\ and 
Z 2 + c ?2 > £ 2 - Clearly, Bl contains B^ We will show how to get a approximation of 
the solution to the PML equation in the truncated domain Bl- 

We introduce local PML problems by using the PML complex coordinate stretch¬ 
ing outside the domain (—ii, Zi) x (C*, Ci+ 2 ). The PML stretching is Xi(x) = (z*, i(*i), 
Xi, 2 (x 2 )) T , which has been proposed in [15] . where Xi t i(x 1 ) = Xi{x\) and 


(2.28) **, 2 ( 3 : 2 ) = 

We define 


*2 + i f** cr 2 (t + Cn+i ~ C,i+ 2 )dt if x 2 > Ci+ 2 , 

*2 if Ci <x 2 < Ci+ 2 , 

X 2 + i J/ 2 CT 2 {t - Ci + Cl )dt if x 2 < ci- 


Ai(x) 


diag 


f x i: 2{x 2 y *i,i(*i) / N \ 

\*i,l(*l) / ’ ii,2(*2)7 ’ 


Mx) = x itl (x 1 yx i> 2(x 2 y. 
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Denote by Df ML = {x = (, X\,X 2 ) £ B L : Q — d 2 < x 2 < (i +2 + ^ 2 }- The local 
PML problems in truncated domains can be defined for some wave source F £ 
_ff 1 (Df ML )' as: find cf> £ 7Tg(f2f ML ) such that 

(2.29) (AiV^jVip) — fc 2 (Jj</), = — {JF,if) V^eEo 1 ^™). 


Then our PSTDDM for the PML equation in a truncated bounded domain consist of 
two steps [3] and |4] 

Before proving the convergence of the method, we introduce some functions and an 
important result which would be used often. The functions are uf, i = 1, • ■ ■ ,N—1, 
and u ~, i = N, ■ ■ ■ , 2 which are defined as: 


(2.30) 

(x) = 

[ + fi+i(y))My) G (xi,yi) d y, * = i,---,iv- 1 , 




(2.31) 

u~ (x) = 

[ ff(y)Ji-i{y)G(xi-i,yi-i)dy, i = N,N~ 1,--- ,3, 

(2.32) 

uf (x) = 

[ iff + fi)Ji{y)G(x 1 ,y 1 )dy. 

J f2iUf22 


The result is about the inf — sup condition of the PML equations in truncated do¬ 
mains. 

Theorem 2.6. Let aod 2 be sufficiently large. There’s some constant a < 1 such 

that 


(2.33) 


sup 

iAeRf(nf ML ) 


(^W, V' 0 ) - fc 2 ( M, ip) ^ .. 0-1/oPMm 

-7T7T7- > T ), 

WWW m inPMLi 1 


where /i _1 < Ck 1+a . C is independent ofk. 

We remark that the recent work (cf. [[T5], 3.16]) of Chen and Xiang shows that 
the inequality in the theorem above holds for a = 1/2.Besides, we know that the 
inf-sup condition number is about k~ l (cf. 132, IT]) for the Helmholtz problem 0 
with Sommerfeld radiation condition (1.2) or Robin boundary condition. 


The proof of the following lemma is omitted. The reader can complete it easily 
by arguments similar to Lemma 3.5-Lemma 3.7 in [ 23 - 
Lemma 2.7. Let 00^2 > 1 be sufficiently large, we have 
(i) For i = 1, ■ ■ ■ , N — 1, 


k + l 


^ Ck(l “)- kL)e 


-■Lik-a 


(ii) For i = N, N — 1, • 


i llif 1 / 2 (dnf“ 1 L ) 


< Ck{l + kL)e-^\ 


iHRBi)' 


1 HRB,)' 


Theorem 2.8. Let <Jod 2 > 1 be sufficiently large, we have 
(i) Fori = 2,--- ,1V- 1, 


f+ - p 

1 1 h- !(f2f ML ) 

(ii) For i = N — 1, • • • ,2, 


< Ck a ^-^k( 1 + kL) 2 e~i k ~< 9 




n-n 




< Ck a{N ~ 


^kfl + kLfe-i^ 
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Algorithm 3 Source Transfer I for Truncated PML problem 

1. Let /+ = /i; 

2. While * * = 1,--- , N - 2, do 

• Find uf € iLg(f]f ML ), where fif ML = {~h -di, l\ +di) x (£* -d 2 , Ci+2 + 
d 2 ), such that 

(2.34) 

(AiVuf, vvo - fc 2 (= (Mft + fi+ 1), </’) W> € ff 0 1 (fi? ML ). 

• Compute T+_ 1 (/+) £ 7L _1 (flf ML ) such that 

^ + +i(/+) = Jr 1 V(A I V(ft+ 1 «+)) + fc 2 (/3+ !«+) . 

• Set /^ 1 = f i+ i + 'F+_ 1 (/+) in fi i+1 n and /+_ x = 0 elsewhere. 

End while 

3. For i = N — 1, find u)y_ 1 £ where = (— 1\ — d\,li + d\) x 

(Cat-i - d 2 ,(N+i + d 2 ), such that \/ijj £ #o( fi jv-i) 

(2.35) (Aat-iVu^^ Vi/)) - k 2 { Jn-iUn-i> VO = (MIn- 1 + /y), V>) • 


Algorithm 4 Source Transfer II for Truncated PML problem 

1. Let /jy = /at; 

2. While i = N, - ■ ■ ,3, 

• Find uf £ 7Lo(fl™ L ) such that 

(2.36) 

(A^iVur, VVO - 

• Compute 'F)li(/, r ) 6 -ff _1 (Il™ L ) such that 

) = tiV^iV^r^r)) + k 2 {fir_ x ur). 

• Set fl~_ 1 = fi— 1 + T“_ 1 (/“) in fVi and f~_ x = 0 elsewhere. 
End while 

3. For i=2, find u 2 £ Ffg(fI™ L ) such that Vip £ Uq(I 2™ l ) 

(2.37) (AiV^, V0) - k 2 (J iu 2 -,^) = (Ji(/ 2 - +/i),^} ■ 


Proof. From the expressions (2.16) and (2.30), we have that uf (x) = uf{x) for 
£ H,; n f2j+i, which implies 


*t + i(f?) = J-'V • (AiV^!«+)) + fc 2 (/3+ ,«+). 
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By simple calculation, we have for any ip G i?Q(fi™ L ), 

lift),*) = -(AiV^,«+ V^)n 1+1 + (AiVu+, V/3+^)n i+1 - </ i+ i,Ptx*) , 

and 

J^+ +1 (/+), V>) = -(A.V/3+ v^)n i+1 + (AiVfit, V/3+ ^)n i+I - fcrf) ■ 

Therefore, 

(j,_ 1 (/+ -;+),«) = - ^ + (/+ r)),«) ninBi 

= -(A i _ 1 V/3 i f (u+ j - ut i), Vu) n . nBi + (Ai-rV^+r - “ti)> V ^ tv )n in ^ 

< Ck 1 ||utr “ ^ti|liji(n 4 nBx.) 

— 11 ^i—i ~ ll_ffi(n? ML ) II v IIh 1 ( 0 F MIj ) ‘ 

On the other hand, uf_ 1 — ut-i = uf_ 1 on <9n™ L and for any ip € #o(f2™ L ) 

(A i _ 1 V(«+ x - utiW) - fcVi-ifot ! ^ uU),ip) = - fti),*) ■ 

By the inf-sup condition |2.33| and Lemma |2.7| we have 


11 ,- 7 + _ || <- r i b 1+a f+ — f+ 

;| u * 1 “i-i ll_f/i(n?_“ L ) — ^ jj-i Ji-i 

+ Cfc 1 + “(1 + kL) IlMi-lUjJl^^QPML) 

<- f+ — f+ 

S ofc /*_! h-l ff - 1(fi p_ T) 




+ C , fc 2+Q (l + /tL) 2 e-2 fc T ff 


Ib^B,)' • 


Therefore, 


ft-Vi 


< Cfc° 


/iti - /+ ! 


_H-i(nf ML ) '* i '* x -ff- 1 (^FA L ) 

+ Cfc 1+a (l + kL) 2 e~^ klS 


IbTb,)' ■ 


(i) follows from the induction argument and the fact that ft ~ ft = 0. Finally, we 
could prove (ii) by an argument similar to that of (i). This completes the proof. □ 
Lemma 2.9. Let 00^2 > 1 be sufficiently large. 

(i) For i = 1,2, ■ ■ ■ ,N — 1, 


ut - ^ + Li(nPML) < Ck m+2 ( 1 + kL) 2 e 


\ 2 f ,-h k 'i< 7 


1 BMBj)' ■ 


(ii) For i = N,-- ■ ,2, 


< Ck a(N - i+ i) + 2(i + kLfe-^ k7(71 


llB 1 (nf“ 1 L ) — V x w & " IIJ llBHBi)' • 

Proof. By using the fact that uf — uf = ut on <9flf ML and for any ip G (fif ML ) 

{Ai V(wt - ut ), VV>) - fc 2 ( Ji(ut ~ ut),ip) = (Ji(ft ~ ft)^) 
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Lemma 2.8 and the inf-sup condition (2.33), the first result can be obtained. The 


second result can be proved similarly. □ 

Theorem 2.10. We define Uq = ufj +1 = 0 in R 2 . Let v = + u~ +1 ) ir 

f 1i D Bl for all i = 1, 2, • • • ,N. Then for sufficiently large &od 2 > 1, we have 

(2.38) [|ii - v\\ h , {Bl) < Ck< N ~^k 2 (l + kLfe~^ 9 \\f\\ m{Bl y , 


where u is the solution to the PML pi'oblem (2.6) in the whole space. 

Proof. From their definitions, we know that uf = uf and u~ = u~ for i = 
1, ■ ■ • ,N. We also define u,} = uT, 1 = 0 in R 2 . Combining with Theorem 


have u(x) = —(u^f 1 (x) + u. 


“JV+l 

■ i+ i(x)) for any x 


€ Lt i 


= !,-■ 


2.5 


we 


, N. Then, by using 


Lemma 2.9 we complete the proof. □ 


Theorem 2.10 shows that the solution v obtained by our PSTDDM is a fine 
approximation of the exact solution to the PML problem in the whole space. We 


remark that the constant ‘C’ in (2.38) generally depends on flf ML and the number N 


of layers due to the inf-sup condition and the induction argument used in the proofs 
which are omitted (cf. [Pd], Theorem 3.7]). 

3. Source transfer block by block. Since every local PML problem in our 
PSTDDM (cf. [l]and[2| is defined on the whole space R 2 , we can use the PSTDDM to 
solve the local PML problems (cf. (2.10) and (2.14)) recursively. As a consequence, 
the domain Bi is divided into some small squares and we only need to solve local PML 
problems defined outside the union of four squares. 

3.1. The PSTDDM for the PML problem in R 2 . We only show the al¬ 
gorithms (cf. [5] and [6| to solve the local PML equation (2.10). In order to show 
the details of our method, Some notations are introduced. For simplicity, we set 
li = I 2 ,, d± = d 2 - In order to use the results obtained in the previous sections, the 
following assumption about the medium property is adopted: 


(3.1) 


/ £-|-cZ/2 nL-fiCL/Z 

01 (t)dt>a, J | 02 (t)dt> 

pl-\-d A-\-d 

/ <Ti (t)dt > <7, / a 2 {t)dt > a. 

Jl+d/2 Jl+d/2 


r l+d/2 


H2 


and 


which is a direct consequence of the assumption HI (2.27). 
We divide the whole space into layers: 


f! 0 = {x = ( X!,x 2 ) T S R 2 : Xi < Ci}, 

LT = {x= {x ll x 2 ) T G R 2 : Ci < xi < Ci+i}, * = 1, • •• ,N, 
Ll N+ 1 ={x= (xi,x 2 ) t G R 2 : Civ+i < xi}, 


and denote by = f2j D LIP 

We define the PML complex coordinate stretching = (x^ J (xi), x 2 J (x 2 )) out¬ 
side the domain (Cj, Cj+ 2 ) x (CiiCi+ 2 ) by x 2 {x 2 ) = x it2 (x 2 ) and 

{ Xl + 1 /c?+2 ai ^ + < ' M+1 ~ 0+2 )dt if xi > Ci+2, 
xi if C j <xi< Ci+2, 

xi + i J^ 1 a lit - Ci + Ci)* if xi < Ci- 


(3.2) x{°{xi) = 
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Then the PML equation’s coefficients are defined as 


A ij(x) = diag 


x\\ Xl y \ 

{x^ix!)'’xY(x 2 y) ’ 


Jij(x) = xY {xi)'xj (x 2 )' ■ 


Let + fi + 1 in fP and = 0 elsewhere. We denote by 7 j + (aq) and 

7 “ ( 27 ) smooth functions such that 7 f(t) = and 7 “(f) = (3~ (t) for any t £ R. 


Algorithm 5 Source Transfer I + for the i th local PML problem 

1. Let fY = fY in fi 1 and = 0 elsewhere, 
while j = 1, • • • , N — 2, do 

• Find u € Lf 1 (K 2 ), such that 

(3.3) -V(A ij -V«+.) - k 2 Jiju+j = JiAf? d + f+ j+1 ), 

• Compute e H _1 (® 2 ) such that 

H j+ i (&) = ^V(A lJ V(7+ +1 n+)) + fc 2 (7+ +1 <,) ■ 

• Set +1 = /+ +1 + 'F+ +1 in fF +1 and /+. +1 = 0 elsewhere. 

End while 

2. Find uf N _ 1 £ 7L 1 (K 2 ), such that 

(3.4) —V(Ai j j\r_iVit) l ) Ar _ 1 ) — k (+ /(yv)- 


Algorithm 6 Source Transfer I for the i th local PML problem 

1. Let f i N = f+ N . 

While j = N, - ■ ■ ,3, 

• Find u~j £ iL 1 (K 2 ), such that 

(3.5) — V(Ai )J -_iVu i j) — k Jij-i'Uij = 

• Compute 'F- ? -_i(/j“-) £ fL _ 1 (R 2 ) such that 

= J yj-i^( A i,j-i^hi-iU~j)) + fc 2 (77-i«ij)- 

• Set in and 1 = 0 elsewhere. 

End while 

2. Find u~ 2 £ Lt 1 (M 2 ) such that 

(3.6) -V(A u Vn- 2 ) - k\J^u- 2 ) = J iA (f~ 2 + fi ji)- 


Algorithm [5] and Algorithm [ 6 ] show the details of our PSTDDM solving the i-th 
PML problem in Algorithm [Tj We omit the details about Algorithm [2] to save the 
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space. Then we can obtain some results similar to those in section [272] but only state 
briefly them when needed. 

The following lemma can be proved by their definitions. We omit the details. 
Lemma 3.1. Let u/f Q = 0 and u~ N+1 = 0. For any x £ &i,j, 

uf{x) = ufj_ i(V) + u~ j+1 (x). 


Lemma 3.2. Let a^d > 1 be sufficiently large. For i = 1, • • • ,7V— 1 

(i) For j = 1, • • • , TV — 2, we have for any x £ ^tj +1 := i x = ( X D X ^) T • Cj+i < 

x\ < Cj +2 and \x 2 - Ci+i| > + d/2}, 

| < Cke-i k ^(\\f\\ H1{Bi y + ||/i|| H i (n . ) ,), 

|Vu+(x)| < Ck 2 e-^(\\f\\ H1{Bi y + 11/^1 * w ). 

(ii) For j = N, - ■ ■ ,3, we have for any x £ 

\^i,j( x )\ < Gke~^ klcr {\\f\\ m(Bi) , + \\fi\\ HHn . y ), 

IV«,~ 3 {x )| < Ck 2 e-^(\\f\\ H1{Bi) , + ||/i|| H1(ni)f ). 

Here 7 = , d =. 

v /d 2 + (2i+2d) 2 

Proof. We give the proof of the first assertion and the second one could be proved 
by the same argument. For j = 1, • • • ,7V — 2, ufj(x) satisfies 

(x) = f f(y)J{y)G{x,y)dy+ f fi(y)J(y)G{x,y)dy 

Jxi<Cj+2 Jxi< Cj + l 

■= ufj(x) + ujf(x) Vx £ Qf j+ 1 - 

It is clear that 

K + ]V)| < C [ \fi(y)\ \G(x,y)\dy + C [ \fi(y)\\G(x,y)\dy, 

7ni». J n»f 

where fl-” = {x = ( X\,X 2 ) T £ : — l — d/2 < X\ < Cj+i} an d Ll/'d = { x = 

{xi,X 2 ) t € f2j : xi < — l — d/ 2 }. By the standard argument (cf. [|T5], 3.11]), we can 
get 

\fi(y)\ \G(x,y)\dy < Ck l/2 e -^ kia \\fi\\ H1{n .y ■ 

i,3 

We recall that Ui(y) and Vwj(y) decay exponentially as \y\ -A 00 , that is 

\ut(y)\<Ck 1/2 e-^\\f\\ HHBl y and | Vu+fo) \ < Ck^e~^ \\f\\ HHBl y 
for \yi\ > Id/2 when ctq d is large enough, which implies 


\fi(y) \\G{x,y)\dy < C sup (+ Vu+ 1 ) 


2 /eo? 


\G{x,y)\dy 


< Che - 1 ** 9 ||/|| ffl 


H 1 (Bi)' ■ 
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Thus for any x € 0^ +1 , 

< Cke-^%\\f\\ m{Bl) , + ||/*L W ). 

It’s so easy to get the estimates, 

\ut]{x)\<Ck l / 2 e-^\\f\\ m{Biy . 

Therefor, we have \ufj(x)\ < Cke~^°(\\f\\ m(Bi) , + \\fi\\ HHn .y) for x £ fI+ j+1 . A 
similar argument implies that 

|Vii+ (*)| < Ck 2 e-^(\\f\\ m{Bl y + H/iL 1 (o.y) € n ti+v 

This completes the proof. □ 

Define In '■= 2l/N+2d and fl™ L := {a; = (xi, X 2 ) '■ Ci~d < x\ < 0+2 +d, (j — d < 
x 2 < Cj +2 + d}. We have the following lemma. 

Lemma 3.3. Assume that a^d > 1 be sufficiently large. There exists a constant 
Cb independent oflN,k and N such that 
(i) For j = 1,2, • • • ,1V — 1, 

II j? i / 2 (an?¥ L ) — Cbk (1 + klN ) ! 


(ii) For j = N, N — 1 


, 2 , 


(3.7) 


\ u i,j II H 1 / 2 (dQ FM \) ~ C' b ^ ^ ^In) ll/ll/fi(_B i )' • 


Proof. The two assertions can be obtained by using arguments similar to those 
in Lemma 3.6 in [15]. We omit the details and show the result 


|u+ 




< Ck 3 ' 2 { 1 + kl N )e-^( ||/||* w + \\fi\\ H1 , ni y). 


By the definition of source transfer operator (2.19), we have 

ii/Hh^bj)' ■ 

Combining the two inequalities above, we have 

< C b k 3 (l + kl N )e~^ | 


|u+.|| 

I llif 1 / 2 ( 8 nf¥ L ) 




Clearly, the C ’s used here are independent of In, k and N. Thus we complete the 
proof of the first assertion and the second one can be proved by the same way. □ 

3.2. The PSTDDM for the truncated PML problem. For the ease of 
presentation, we denote 


ffi u = (Ci-rf,C 2 ) x (0-00+2 + d), 

= ( 0 , 0 +i) x (o - d, 0+2 + d), j = 2, • • •, iv - 1, 
U= (Cn , Ov+i + d) x (0 — d , 0+2 + d). 
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We can get the approximation uf(x) of uf(x), i = 1,- •• , N — 1, in flf ML by 
using Algorithm [ 7 ] and Algorithm [s] . where fit(ar) are defined by 

(3.8) u+ = fi+ in fit = u~ N in 

uf = + K,j+i in for j = 2, • • • , N - 1. 

In Algorithm [7] and Algorithm |8] we have defined ft,j = ft + /*+1 in n Tj and f+ 
are defined by 

1. Let /+ = /1 in fl™ L . 

2. Compute G dd _1 (nf ML ) such that 

^+ +1 = Jr 1 V(A i V(/3+ ,«+)) + k 2 V3+ +1 &t) . 

3. Set ff + 1 = / i+ i + in fi i+ i n B L and ff +1 = 0 elsewhere. 

We also could obtain the approximation u~ (x) of u~(x), i = N, ■ ■ ■ ,2, in fl™ L . The 
details are omitted in order to save space. 


Algorithm 7 Source Transfer I + for local Truncated PML problem 

1. Let /+ = f+. 

While j = 1, • • • , N — 2, do 

• Find fit. g #0 where = (C i~d, Ci +2 +d) x (Q-d, 0 + 2 +d), 

such that Vi/> € tdo(^™ L ) 

(3.9) (AjVfit Vi/>) - ^(tijfit.,^) = (Ji(/+ + /t. +1 ),V>) , 


• Compute ^t. +1 (/+.) g H 1 (flP“ L ) such that 


if/+ 


(/+) = , j V (A j i j V ( 7 1|., fi t^.)) + fc 2 ( 7 t + 1 fi+.) 


• Set f+ j+1 = /t . +1 + T+- +1 (/+) in D j+1 n nff L and f+ +1 = 0 else¬ 
where. 

End while 

2 . Find u+ N _ 1 e where = (C* - d, Ci +2 + d) x ( 6 v-i - 

d, Cy+i + d), such that Vi/’ G H/J (^™-i) 


(3.10) 

(Ai ! AT_lVfit JV _ 1 , V-0) — fc = ^dj i iV-l(/j|jv_ 1 t ■ 


We can improve the local inf-sup condition (2.33). 

Lemma 3.4. For sufficiently large a$d > 1, we have the inf-sup condition for any 

<t> G ddo(^™ L ) 


(3.13) 


(A iyJ V<f,Vip) - k 2 (J itj f>,ip) ^ IUII 
SU P -jmij- >MlMl ff i(QeM L) , 


where p 1 < Ci S (l]yk ) 3 ^ 2 if lisik large enough, and p 1 < Cis i/ /at/c « 1 . Ci S 
independent ofli\r,k and N. 
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Algorithm 8 Source Transfer I for local Truncated PML problem 

1. Let f i M = ffa. 

While j = N, - ■ ■ ,3, 

• Find G suc ^ that \/i/j G 

(3.11) (Aij—i\7u i j,S7'ip) — k , 


• Compute 'F ij _ 1 (/j j ) G H such that 

fa-i (/■■) = +k 2 ( 1 T_ 1 u7 j ). 

• Set fr.^ = /+_ 1 + vE>~ / _ 1 (/~ / ) in and = 0 elsewhere. 

End while 

2. Find «r 2 g fa(flffa) such that G fl™ L ) 

(3.12) (A u Vu- 2 , VVO - fc 2 (J u «- 2 ,^) = (Ji,i(^7 2 + ftM) • 


Proof. The inequality can be proved easily by using scaling argument. We know 
that there is a unique solution <f> G to the problem 

-V(AjjV^) - k 2 Jij4> = F, 

for some F G H^Alffa) 1 . We define a mapping m : I := [0,1] x [0,1] —> as 

m(z) = Inz + (Ci — d,fj ~ d) and denote by <f>(z) := <j>(m(z)) and F(z) := F(m(z)). 
The equation above implies <j>{z) G Hfa) satisfying 

(3.14) -V z (A itj (m(z))V z fa)) - ( l N k) 2 Jij(m(z))fa) = l 2 N F(z). 

If Ijsik large enough, by the local inf-sup condition (2.33[), we get 


(3.15) (( l N k ) 2 fa) 


LAP 


fa 


m(i) 


where |H|#i(/)/ is defined as 


sup 


) 1/2 <C is {l N kf' 2 l 2 N F{z) 


fa 


HHiy 


^HAI) (( l N k ) 2 ||V»Hi2 (/) + 


1/2 ’ 


from the definition of weighted norm ||u|| H i( u) at t he beginning of section 
ever, if l^k ~ 1 it’s known that the problem (|3. 14) is elliptic, then we have 


2.3 


How- 


(3.16) 


( fa) 


lai) 


fa) 


) 1/2 <C ls l 2 N F(z) 


hai) 


HAiy 


Clearly, C zs is independent of In, k and N. Finally, The consequence is obtained by 
combining the inequalities (3.15), ( |3.16| ) and the fact that 

„ 2 2 „ 2 2 

(l N k) 2 fa) __ = k 2 ||<^(a;)|| i2(n pML ) , fa) = Ifa^m^™^ , 


l 2 N F{z) 


hai) 


LAI) 














































18 


Y. Du & H. Wu 


□ 

In general, we can expect that l^k is less than k. If N is large enough such that 
lisrk ~ 1, the local truncated PML problems (cf. (3.9) (3.12)) needed to be solved are 
about elliptic. 

Lemma 3.5. Let a^d > 1 be sufficiently large. There are constants Ct and C b t 
independent of In , fc and N such that 
(i) For i = 1,2,- •• ,1V — 1, 


(3.17) 


lift — itt I 


"i llj? 1 (flf ML ) 
(ii) For i = N,N — 1,- ■ ■ ,2, 
(3.18) 

||^i _ Ui ||ifl(fiPML 

Here Ck,N,j, j £ N are defined as 


<C bt C KN>i k 3 {l + kl N ) 2 e-i k ^\ 


< CbtCk,N,N+i-ik 3 (l + fc/jv) 2 e 2 fc7IT 


(3.19) 


j / -V —i 

cw = E(£ ( c ^~ nP 

q=l \ 






Proof. We only show the details of the proof for the first assertion and the second 
one could be proved similarly. At the beginning, we recall the property (cf. |15j . 
Theorem 3.7) of source transfer operators that there’s a constant Ct independent of 
In, k and N, such that 


f ~ f 

J l,J J l,J 


H 1 ( QPML), - Ct \\ U ij -1 ’ 


L-7+ 


\fi fi llffi(fif ML )' — Ct z b-i||i/i(n i nB t ) ’ 


for i, j = 2, • ■ ■ ,N — 1 from their definitions and calculations similar to (2.19)-(2.20). 


Using the argument in Lemma 2.9 and Lemma [3.3| , it’s easy to get 

(3.20) ||G+. - fit.11 < p~ l (f+ + f+ +1 ) - (/+ + /+ +1 ) 




T P 1 (1 + kljy) 11 ut 






II f+ 


ft-ft 


lffi(n™ L )' 


+ C b p~ l k 3 {\ + kl N fe-^-° \\f\\ m{Bl y . 
By the induction argument and the fact that 


%ill 


<^ _1 n ft-n 




+ C b p- 1 k 3 (l + kl N ) 2 e~^ ff \ 


I H^(B t y > 


(3.20) implies for j = 1, ■ ■ • ,N — 1 


l-i 


(3.21) 


\ U i,j ^*,jllfl’ 1 (fif¥ L ) 


<J2(. C tT~ 1 ) P -T- 1 - 


p =o 


[ || ft- /t Ll (n PMLy + C b k 3 ( 1 + kl N fe ||/||* w ] 
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Similarly, we have for j = N,- ■ ■ ,2 


(3.22) 


I U i,j 




I ft-ft 


i \\H 1 (Qf M ^) 


p—0 

,+C b k 3 (l + kl N ) 2 e~^ 9 \ 


From (3.21), (3.22), Lemma 3.1 and the definition 3.8 we obtain 

N—2 


(3.23) 


I u t ^llj?i(nf ML ) 


■ 1 '' p ■ fi 1 - 


< E ( c * 

p—0 

ft - ftWm^y + C ^0 + kl N ) 2 e~^ 9 \ 
< Y, (Q/x- 1 ) p [||'u+i-ut rll^^nB^ 


N—l 


P =1 


+ ^k 3 {l + kl N fe-^\\f\\ m{Bl y]. 

Since C t and C b don’t depend on lw,k and N, we can denote C bt = 0). Then we 
complete the proof for the first assertion (3.17) by the induction argument and the 
fact 

N—l 

K - fif || ffl(n PML } < E + kl N fe\\f\\ m{Bl y . 

P=1 


\w(B t y J • 


J • 


□ 

The following theorem is a direct consequence of Theorem |2.5[ Lemma |3.5| and 
the fact that uf 1 = uf in fl, ; U fb+i. 

Theorem 3.6. Letu q" = f% + i = 0 m and u(x) = —(iif_ 1 +u~ +1 ) infliHBL 
for all i = 1, • • • ,N. Denote Ck t N = Ck t N,N-i- Then for sufficiently large a^d > 1, 
we have 

(3.24) ||u — '^f\H 1 {B L ) — C b tCk,Nk 3 (l + klN) 2 e 2 fc 'i' cr ||/|| ff i( Sj y . 


We remark that from the theorem above, we can know that the larger number N 
doesn’t mean the solution u performing better although the local problem solved may 
be elliptic. However, our numerical examples in the following section show that the 
relative errors between u and the discrete solutions don’t increase significantly when 
N becomes larger. 


for constant and heterogeneous wave number by FEM and STDDM, where / is given 
so that the exact solution is 

— r 3 (r 3 + 3r 2 — 12r + 9)Hjf\kr), r < 1, 

— H^\kr), r >= 1. 

Let di = 0.2, d 2 = 0.1 and h = I 2 = 1.1. So the computational domain B b is 
(—1.3,1.3) x (—1.2,1.2) and the perfect matched lay is Bjf\Bi where Bi = (—1.1, l.l) 2 . 
It is easy to see that u £ C 2 (R 2 ) and supp / c B;. 



4. Numerical examples. In this section, we simulate the problem o and|1.2 
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We define the medium property m by setting l\ = l 2 = 1.18 and (Jj{t) = 
aj(t) + (t — lj)a'j{t) for lj < t < lj, where 


(4.1) 


&j(t) = 7o 





h)\h 



The functions (3f(x 2 ), x 2 6 Hj, i = 2, ■ ■ ■ ,N — 1, used in the source transfer 
algorithm are defined as 


( 1, Ci < x 2 < C i + Ci + AC/4, 

^(* 2 ) = < r)i{x 2 ), 0 + AC/4 < a: 2 < Ci + 3AC/4, 
l 0, Ci + 3A//4 < x 2 < Ci+i, 


and /3j = 1 — /3/~, where 


??i(z 2 ) = 1 + 


( X 2 - (Ci + AC/4) 

V ac/2 


/ x 2 - (Ci + AC/4) 

l AC/2 


2 


Clearly, (a; 2 ), i = 2, - ■ ■ ,N — 1, are in C 1 (fij) and this fact avoids the discontinuity 
of j3f(x 2 )' which may make ff~ oscillate heavily. 

We use the finite element method to solve truncated PML problems. The number 
of nodes in the ary-direction is rij = q ■ 2Lj/\, j = 1,2, where q is the mesh density 
which is the number of nodes in each wavelength A = 2ir/k. Then the number of 
degree freedom DOF is nin 2 . Let N be the division number in the a; 2 -direction. e*, 
ef and e s denote the relative error in ffCseminorm of the interpolation, the FEM 
solution and the PSTDDM solution bounded in Bi respectively. 

We first test the algorith m [3|a nd|4|for the wave number k/(2n) = 25 and k/(2n) = 
50. The left graph of Figure |4~l] plots the relative error decay of the interpolation, 
FE solution and PSTDDM solution with a fixed number of lays N 2 = 10 in terms 
of DOF for k/(2n) = 25, 50 respectively. We could find that the relative errors 
of PSTDDM solution is the same to that of FE solution when DOF is equal. This 
is best result about comparison between the PSTDDM and FEM which we could 
expect, since the details of the algorithms [3] and [4] show that the errors of PSTDDM 
solutions can not be less than those of F E solutions under the condtion that the mesh 
is same. In the right graph of Figure [ tT] we set DOF = 624 x 10 4 and give the relative 
errors in H l -seminorm of the PSTDDM [3J [^solutions in terms of the number of lays 
in a; 2 -direction N = 1,5,10,20,25,50,100, for k/(2n) = 25,50 respectively, where 
N = 1 means that this solution is the FE solution. It is shown that the error of 
PSTDDM solution remains unchanged even if the number of lays in the ® 2 -direction 
becomes larger. So we could choose a relatively large number of lays to reduce the 
computational complexity. 

Next we test our further consideration (cf. |3.8| [7J [8]) about the PSTDDM for 
k/(2ir) = 25 and k/(2n) = 50. The parameters about PML layers are still those 
provided at the beginning of this section, since they’re not essential from the previous 
proofs. 

In the left graph of Figure |4.2| we set N = 10, and show the error decay of 
the FE solution and further PSTDDM solution when mesh density q increases. The 
graph is very quite similar to that of Figure |4~T| what we could like to obtain. In the 
right graph, we show the relative errors of the further PSTDDM (cf. |3.8| [?J [8]) when 
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N= 5,10, 20, 25. Thus the number of the squares, which the domain Bi is divided into, 
is N 2 = 25,100,400,625. 

We remark that it’s not necessary to set the number of layers N too small because 
of the fixed width of PML layer resulting in low computational efficiency in practical 
application. 




Number of lays 


Fig. 4.1. Left graph: The relative errors e^, ef, e s for the interpolations, FE solutions and 
PSTDDM [3] 2| solutions with a fixed number of lays in X 2 ~direction N = 10 in terms of the number 
of degree freedom DOF = niri 2 for k/(2n) = 25 and k/(2n) = 50 respectively. Right graph: The 
relative errors for the PSTDDM [3j 2| solutions in term of the number of lays in the X 2 ~direction for 
k/( 27r) = 25 and k/(2 tt) = 50 respectively, and setting DOF = 624 x 10 4 . 




Fig. 4.2. Left graph: The relative errors ef, e s for the FE solutions and further PSTDDM 
(cf. 0 solutions with a fixed number of lays N = 10 in terms of the number of degree freedom 
DOF = nin .2 for fc/(27r) = 25 and k/(2n) = 50 respectively. Right graph: The relative errors for 
the further PSTDDM (cf. \3.8\ [?j [3j )p solutions in term of the number of squares for k/(2n) = 25 
and k/{2 , K) = 50 respectively, and setting DOF = 624 x 10 4 . 
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